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Abstract 

At HERA, the electron-proton collider at DESY, Hamburg, the large flux of 
almost on-shell photons accompanying the lepton beam is being used to shed new 
light on the structure of the photon. Recent results are reviewed and discussed, with 
emphasis on those aspects of the photon's nature which should be understandable 
using perturbative QCD. 
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1 Introduction 



At DESY, Hamburg, 820 GeV protons collide with 27 GeV electrons (or positrons) in 
HERA. Two general purpose detectors, HI [ffj and ZEUS 0, are positioned at opposing 
interaction points on the accelerator ring. The high flux of almost on-shell photons which 
accompanies the lepton beam provides a unique opportunity to study the nature of the 
photon and its interactions. 

Collisions take place between protons and almost on-shell photons (with a negative mass 
squared in the range lCT 10 < Q 2 < 4 GeV 2 ), at •yp centre-of-mass energies (W JP ) up to 
300 GeV, an order of magnitude greater than previously available. HERA does not, of 
course, measure the structure function F%, as this is defined in terms of the scattering 
of electrons off photons. However, studies of jet production in yp interactions provide 
information about the photon structure down to momentum fractions as low as 5 x 10 -2 , 
and because the photon is being probed by partons from the proton, there is a high 
sensitivity to the distribution of gluons in the photon. 

The presence of a 'hard' energy scale (usually from the high transverse energy (E^ 1 ) of 
the jets) means that perturbative QCD calculations of event properties can be confronted 
with experiment. Examples of leading order (LO) QCD diagrams for photoproduction 
at HERA are shown in figure [l]. At this order, two processes are responsible for jet 
production. The photon may interact directly with a parton in the proton (figure |T]a), 
or it may first fluctuate into an hadronic state (figure [l|b). In the first case, known as 
the direct contribution, the fraction of the photon momentum (x 7 ) participating in the 
hard process is one. The final state consists of two jets, the proton remnant and the 
scattered electron. In the second case, known as the resolved contribution, the photon 
acts as a source of partons, which then scatter off partons in the proton, and x 1 is less 
than one. The final state in this case includes a photon remnant, continuing in the 
original photon direction, in addition to two jets, the proton remnant and the scattered 
electron. Direct photon events do not involve the photon structure, but can probe the 
parton distributions in the proton down to x p ~ 10~ 3 , where x p is the fraction of the 
proton's momentum entering into the hard process. This process is directly sensitive to 
the gluon distribution in the proton, and complements indirect extractions || using the 
measurement of F 2 in deep inelastic scattering (DIS). 

Seen from a perspective of understanding and testing QCD, photoproduction at HERA 
provides information about quark and gluon interactions in the yp system over a broad 
range of kinematic variables, information that is complementary to that obtained from 
DIS, either in ep events at HERA or in ej events at TRISTAN and LEP. 

Both the ZEUS and HI detectors have high resolution calorimeters and tracking chambers. 
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Jets have been measured in the pseudorapidityQ region — 1 > r\ > 2.5 and with E 3 ^ 1 in the 
range 6 < E J T et < 37 GeV. 





a) b) 



Figure 1: Examples of leading order diagrams for a) direct and b) resolved photopro- 
duction. 



2 Untangling the Subprocesses 



Because of the large cross sections involved, photoproduction physics provided interesting 
results very soon after HERA began colliding the beams. At that time, evidence for 
'hard' interactions in the resolved photon process had just been observed for the first 
time, by the AMY collaboration Both the ZEUS and HI collaborations rapidly 
provided confirmation of hard scattering from partons inside the resolved photon, this 
time in 7^ events [||]. Jet production was seen, and energy flow in the rear (electron) 
direction provided clear evidence for the presence of a photon remnant. 

In events with two (or more) jets, the parton kinematics of photoproduction events can 
be reconstructed. In particular, the momentum fractions of the photon and proton can 
be measured. For two-to-two parton scattering in LO QCD, energy and momentum 
conservation give the fraction of the photon energy involved in the hard scatter as 



Erpparton _ v parton 
partons f^T 6 /-,\ 

7 " 2yE, 



1 The z axis is defined to lie along the proton direction, and rj = — ln(tan |) where 9 is the angle 
between the jet and the z axis. 
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partons e /q\ 

Xp ~ 2E P ' [) 

where yE e is the initial photon energy, E p is the proton energy, and the sum is over the 
two final state partons. For direct photon events, x^° = 1. For those events in which 
two jets are observed in the detector, it is possible to estimate the fraction of the photon 
momentum which participates in the hard process by replacing the sum over partons in 
equation ([!]) with a sum over jets. It is then possible to make an explicit definition of the 
observable x~ 
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where now the sum runs over the two jets of highest E^ 1 . An analogous definition can 
be made for the proton momentum fraction. In the x® BS distribution thus obtained, 
the LO direct and resolved processes populate different regions, with the direct processes 
concentrated at high values of x° BS . The peak arising from the direct contribution will 
not necessarily lie exactly at x° BS = 1, due to higher order effects and/or hadronisation, 
but will still correspond to the kinematic region where most or all of the energy of the 
photon participates in the hard subprocess. This treatment avoids possible confusion 
arising due to the fact the simple separation between resolved and direct processes is 
only possible at LO in QCD and x^° becomes ambiguous even at next-to-leading order 
(NLO). 

The x® BS and x p BS distributions (not corrected for detector effects, and therfore called 
x meas ) measured by the ZEUS collaboration using 1992 data are shown in figure |2|. 
The x° BS distribution (figure ^b) rises at both low and high values. The Monte Carlo 
simulations of the resolved and direct processes (which also include detector effects) shown 
in the same figure, have very different characteristics. The resolved processes show a rise 
towards low x° BS , as observed in the data, but cannot account for the rise at high x° BS . 
The direct processes predict a sharp rise towards high observed in the data and 

only a small number of events for x® BS < 0.7. The conclusion is that the peak at the 
high end of the x® B distribution results from direct processes. 



3 Jet Cross Sections 



The results presented in the previous section show good qualitative agreement between 
the data and QCD expectations, in the sense that the expected production mechanisms 
and event topologies are observed. However, signs of discrepancies are already seen. In 
particular, when the transverse energy flow around the jets is studied, more energy was 
seen in the data in the forward (proton) 
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Figure 2: Kinematic distributions for events with two or more jets, (a) x™ 
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region than was expected from Monte Carlo simulations J/]. In addition, in figure 
there is an excess of low x° BS events in the data over the simulation. In order to obtain 
quantitative information from photoproduction processes, clearly defined cross sections 
must be measured and compared with theoretical calculations. I concentrate here on jet 
cross sections, although much information can be also be obtained from measurements of 
charged particle spectra Hj. 



3.1 Inclusive Jet Cross Sections 

Inclusive jet cross sections provide a good test of perturbative QCD, and have been 
measured by both the HI |?j and ZEUS collaborations. At both experiments, jet 
finding is performed using a cone algorithm The cone radius of the jet finding 



algorithm is typically defined such that R = J Ar] 2 ell + A0 2 eZ/ = 1, although other values 
could also be used in the future. The algorithm searches in pseudorapidity-azimuth [r\ ce u- 
4>ceii) space for the cone containing the highest summed Et- 

HI use the electron calorimeter of their luminosity monitor to measure the energy of the 
low angle scattered electron in the range 0.25 < y < 0.7, where y = E 1 /E e , the fraction 
of the initial electron energy (E e ) carried by the almost real photon with energy E 1 . The 
corresponding range of photon virtualities is approximately 10~ 8 GeV 2 < Q 2 < 0.01. 
For jet studies at ZEUS, photoproduction events are usually defined by demanding that 
the electron is scattered at small angles and does not emerge from the beam pipe with 
the central ZEUS detector. This requirement corresponds approximately to a cut of 
Q 2 < 4 GeV 2 . The range of y values used by ZEUS is 0.2 < y < 0.85. 

The inclusive jet cross cross section da/dt] measured by the ZEUS collaboration is shown 
in figure ^. Also shown are the cross sections calculated by the PYTHIA [|T1|] Monte 



Carlo model (including the effects of hadronisation and parton showering). The model 
with both the photon parton distribution sets shown describes the data well, with only 
a small discrepancy in the foward (high r/) region for the low E^ cross section. The 
inclusive jet cross section da/dE^ 1 is shown in figure |], for two different ranges of r?- 76 *, 
one of which includes the forward region and the other of which excludes it. 

Given the failure of the Monte Carlo models used here to describe the energy flow in 
the forward region shown in ref. 0) it is not surprising that the cross sections do not 
completely agree in the forward region at low E^ 1 . At higher the uncertainties 
arising from hadronisation can be expected to become less significant, and any test of 
QCD via such comparisons becomes less ambiguous. NLO calculations of these processes 



are already becoming available [12 . 
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Figure 3: Differential ep cross section da/drj for inclusive jet production, integrated 
over Ej? from three different thresholds (E^ > 8, 11 and 17 GeV) in the kinematic 
region defined by Q 2 < 4 GeV 2 and 0.2 > y > 0.85. PYTHIA calculations are shown for 
comparison. 
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Figure 4: Differential ep cross section da/dE^ 1 for inclusive jet production, integrated 
over two rj^ et ranges (— 1 < rf et < 2 and — 1 < r/- 76 * < 1) in the kinematic region defined 
by Q 2 < 4 GeV 2 and 0.2 > y > 0.85. PYTHIA calculations are shown for comparison. 
For the latter 77 je * range, both Monte Carlo curves and data have been multiplied by 0.1. 
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3.2 Dijet Cross Sections 



Although inclusive jet cross sections have the advantage of calculational simplicity, there 
are advantages in measuring dijet cross sections. In this case, reconstruction of the kine- 
matics of the hard scatter becomes possible, in particular of x° BS and x® . By making 
kinematic cuts it is possible to untangle some of the many unknowns in the photoproduc- 
tion process, in particular the dependence upon the parton distributions in the proton 
and photon. 

The ZEUS collaboration has measured the differential cross section da/dfj for dijet pho- 
toproduction Jl3| 

(ep — ► e7 p —> e X + 2 jets) in the y interval 0.2 < y < 0.8, for events with at 
least two jets of E^ 1 > 6 GeV, and a difference in the pseudorapidities of the two jets 
of |A?7| < 0.5. The variable fj = |(r/i + r] 2 ) is the average pseudorapidity of the two 
jets of highest transverse energy. This cross section is measured separately for resolved 
and direct photon events, as defined by a cut on x° BS at 0.75, with the low x° BS events 
identified with resolved photoproduction and the high x® BS events identified as direct 
photoproduction. The expression for x® BS (equation [|) can be rewritten in terms of |Ar;| 
and fj, assuming the jets to have equal transverse energy, to obtain 



A similar expression can be written for the proton, 



^ = ^cosh^. (5) 

Constraining the jets to be at almost equal pseudorapidities, with the cuts on |A^|, means 
that the hyperbolic cosine term takes its minimum value of unity. Thus the minimum 
available x values are probed for a given E^ and fj, and there is a strong correlation 
between fj and x° BS in the direct cross section and between fj and yx® B in the resolved 



cross section [ 14] . 



The direct cross section is shown in figure [5]a, compared to leading order QCD calculations 
using various different parton distribution sets for the proton. The parton distribution 
sets used are the GRV LO |5|, CTEQ2M |B|, MRSA |T7| and MRSDq parton 
distribution sets for the proton. The shape of the direct cross section differs from that 
of the LO QCD calculations. However, several effects (apart from the choice of parton 
distribution set) influence the comparison between data and theory. QCD calculations 
are only available at LO, and higher order corrections may be large. Also, due to the fact 
that we are probing low x partons in the proton, the standard approximation that the 
incoming partons are collinear and on-shell may be invalid |2Uj . Finally, non-perturbative 
'hadronisation' effects can be expected to be significant. Considering these uncertainties 
the agreement with the theory is fairly good. 



8 




Figure 5: da/dfj for ep -> eX + 2 jets, | Ar?| < 0.5, £^ e * > 6 GeV, 0.2 < y < 0.8, Q 2 < 
4 GeV 2 , a) x° BS > 0.75, b) x° BS > 0.75, The solid circles are corrected ZEUS data. The 
inner error bars indicate the statistical errors, the outer error bars show the systematic 
uncertainty excluding the the correlated uncertainty from measurement of energy in the 
calorimeter and the integrated luminosity, which are show as a shaded band. In a) the data 
are compared to LO QCD calculations using several parton distribution sets for the proton 
and the GS2 set for the photon. In b) the data are compared to LO QCD calculations 
using the LAC3, GS2, GRV, LAC1 and DG parton distribution for the photon. The 
proton parton distribution set used is the MRS A set. 
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Figure [5|b shows the measured cross section da/dfj for resolved photoproduction. The 
LO cross sections, shown for comparison, are calculated using different photon parton 
distribution sets and the MRSA [[17] set for the proton. The theoretical sensitivity to 
the parton distributions in the proton is small (not shown), with the variations between 
curves calculated using different parton distribution sets being much less than the esti- 
mated errors on the measured cross section. The DG pif , GRV [TS| and GS2 |f22[ parton 
distribution sets reproduce the shape of the cross section well and can be brought into 
agreement with the data by applying a multiplicative factor of around 2. Factors of this 
size may reasonably be expected to come from NLO calculations []12|. The LO calcu- 



lations using the LAC1 and LAC3 |23j parton distribution sets cannot be brought into 
agreement with the data by a constant normalization factor. 



4 Photon Structure 

The Monte Carlo simulations discussed so far attempt to describe photoproduction by 
a hard parton-parton subprocess coupled with parton distributions in the proton and 
photon and parton shower and hadronisation models. No allowance is made for the 
possibility of further interactions (other than a single parton-parton scattering) taking 
place in a 7p event. In obtaining the results shown, the simulations are only used to 
correct for detector acceptance and smearing. The effects of the discrepancies between 
the Monte Carlo simulations and the data on the measurement of the cross sections are 
estimated and included in the systematic errors. Thus the cross sections can be presented 
as measurements of observed jet cross sections with confidence. However, this confidence 
evaporates when we come to interpret the cross sections in terms of parton distributions. 
If there are unknown other processes contributing to the energy flow in the event, they will 
distort the correlation between 'parton' kinematic variables in the model and jet variables 
in the data. 

As an example of this, figure [| shows the uncorrected transverse energy flow per jet 
l/NdEx/ d5rj around the jet axis for events entering the two dijet cross sections of fig- 
ure |. For the events with x^ 1 > 0.75 (direct events, figure [|a and figure ||a), both the 



HERWIG [33] and PYTHIA simulations reproduce the data distribution well. The same 
distribution is shown in figure [|b for events with x^ al < 0.75 (resolved). In this case both 
simulations fail to describe the pedestal in the energy flow in the forward region, as was 
also observed (for inclusive jets without an x° BS cut) in references and ||. Thus when 
interpreting the cross section in figure |5|b in terms of parton distributions, it is difficult to 
estimate what the effects of hadronisation and parton showering might be on the shapes 
of the curves, since the jet pedestal energies are not well described in the Monte Carlo. 
Therefore no parton distribution set (with the exception of LAC3) can be confidently 
excluded at this stage. 
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ZEUS 1993 




Figure 6: Figures a) and b) show the uncorrected transvserse energy flow dE T /d8r] 
around the jet axis, for cells within one radian in <fi of the jet axis, for a) direct and 
b) resolved events. The solid (dashed) line represents the distribution from PYTHIA 
(HERWIG). 
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4.1 Multiple Interactions 



One candidate for improving the agreement between the Monte Carlo models and the data 
is to allow for the possiblity of more than one hard interaction per 'yp event [25|. This 
can happen because the parton densities probed at HERA are sufficiently large that even 
from a naive probabilistic point of view, the chance of more than one hard scatter in a 7^ 
interaction becomes significant. Models of this nature are usually called 'eikonal' models 
and are also used to regularise the strong growth exhibited in minijet calcuations of the 
7P cross section 
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It has been demonstrated by HI [|27] that including some form of 
multiple interaction model can bring the Monte Carlo simulation into significantly better 
agreement with the data. As multiple interactions in this type of scenario are generated 
by the density of partons in the photon, they are expected to be absent in the direct 
process, where the photon parton distributions play no part. Returning to the jet profiles 
in the resolved and direct samples isolated by ZEUS in reference |L3[ (figure the fact 
that the forward energy excess is not observed in the direct (high x® B ) events lends 
support to the idea that multiple hard interactions are a reality. 



4.2 Parton Distributions 

It is possible to use the PYTHIA multiple interaction model, together with an additional 
'underlying' jet pedestal energy and the usual parton showering and hadronisation models, 
to obtain a good description of the energy flow around the jets jP7| . Within this model, 



there then exists a particular correlation between the LO, generated partons and the 
observable jets. In the absence (so far) of better theoretical tools, it is interesting to 
assume that the Monte Carlo model represents the 'truth' as far as this correlation is 
concerned, and deduce what parton distribution is needed in the Monte Carlo simulation 
in order to successfully describe the observed jet rates. This procedure was carried out 
by HI, and the resulting distribution is shown in figure [7[ 

Once this has been done, the fact that measurements of constrain the quark distri- 
butions in the photon in this kinematic region can be used to deduce what part of the 
parton density is attributable to gluons in this model. The result is shown in figure M. 



4.3 Remnant Structure 

Most analyses of photoproduction events at HERA have implemented a cone algorithm to 
find jets. The cone algorithm uses a cone of fixed radius in pseudorapidity and azimuthal 
angle space and is well suited for high transverse energy jets at fairly central pseudora- 
pidities. However, the photon remnant will mainly deposit energy at low angles relative 
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to the electron direction, over a large negative pseudorapidity range, and a clustering 
algorithm is more appropriate. ZEUS have used the kx jet clustering algorithm in the 



laboratory frame |28|, |29[ to isolate and study this photon remnant |30 . 

This algorithm finds jets by calculating a quantity for each calorimeter cell 'hit' with 
respect to all other hit cells. In the small angle approximation, kx reduces to the transverse 
momentum squared of the lower energy cell with respect to the higher energy cell. The 
value of kr for each pair of cells, or clusters, is calculated from the formula: 

k T = 2mm(E 2 n , E 2 J(1 - cos# nm ), (6) 

where E n and E m are the energies for cells n and m, and 9 nm is the angle between the two 
cells. When all of the kx values of an event have been calculated, the two cells with the 
lowest kx value are merged. This process is iterated, merging cells or clusters of cells, until 
the desired number of clusters remain. The proton remnant is included in the clustering 
procedure as a pseudo-particle which moves in the proton direction. For high transverse 
energy jets the K t cluster finding algorithm gives results similar to those obtained with 
the usual cone algorithm. 

The kx algorithm was forced to find three energy clusters in each event, in addition to the 
cluster related to the proton remnant. The pseudorapidity r] clusteT distributions of the three 
clusters obtained with the kr algorithm are shown in figure |S]a-c, in order of decreasing 
transverse energy. The data ( full circles) and MC ( resolved and direct, histogram) are 
shown normalized in the region rj cluster < 1.6. While the first and second clusters are 
predominantly found in the rf luster > region, the third cluster has an accumulation in 
^duster < ^ e the outgoing electron direction. The agreement between the data and 
the MC is reasonable for all three clusters except for the very forward region r) cluster > 1.6 
where the discrepancies already noted in the forward energy flow can be expected to be 
significant. In particular the excess observed in figure |c in the negative pseudorapidity 
region is well described by the Monte Carlo simulation including resolved and direct 
processes (full histogram). However the (LO) direct process does not contribute to this 
excess as shown by the dashed line in the figure; its distribution is rather flat as expected 
from the absence of photon remnants for this class of events. As a consequence, the third 
cluster in the negative pseudorapidity region can be identified with the photon remnant. 

After requiring the first two clusters to be jets with E T et > 5 GeV and pseudorapidity 
rf et < 1.6, as in reference ||, the rj distribution of the third cluster is shown in figure 
Id. Compared with the MC expectations, the data show an excess toward higher rj 
values corresponding to higher angles with respect to the electron direction. A cut of 
^clusters < _y was then applied to the third cluster. After these cuts, x® BS was found to 
peak at low values as expected for resolved processes. 



The energy distribution for the third cluster is shown in figure [T0| a The solid histogram 
shows the Monte Carlo expectation which agrees reasonably with the data. The trans- 
verse momentum (with respect to the beam) distribution of the third cluster, shown in 
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figure |lO]b, peaks at 1 GeV with a tail extending to 6 GeV. The Monte Carlo expectation, 
shown in histogram, shows a lower average value. This is in qualitative agreement with 
the expectation is that the 'anomalous' component of resolved photoproduction, where 
the photon splits into a perturbatively calculable qq pair, should have a higher trans- 
verse energy remnant than the so-called Vector-meson (VMD) like component, where the 
splitting is non-perturbative. The Monte Carlo models so far used assume the VMD-like 
topology for all events. 



5 Summary 

HERA data provides unique opportunities to study the structure and interactions of the 
photon. Jet production by both the whole photon and its parton constituents have been 
unambiguously observed, and jet cross sections are being used to test QCD and constrain 
the parton (particularly the gluon) distributions in the photon and the proton. Studies 
of the structure of the photon remnant give indications of a higher transverse remnant 
energy than is expected in simulations which describe the photon as a vector meson. 
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Figure 7: The distribution unfolded using PYTHIA from events with two jets 
E J T et > 7 GeV and if et < 2.5. 
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Figure 8: LO estimate of the gluon density in the photon unfolded using PYTHIA from 
events with two jets E 3 T et > 7 GeV and rf et < 2.5. 
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Figure 9: Pseudorapidity distributions of the first (a), second (b) and third (c) clusters 
found by the Kt algorithm: ZEUS data (full circles), MC direct +resolved (full histogram) 
and MC direct (dashed histogram), (d) as (c) but imposing a cut at P\ 



clusteri.2 
T 



> 5.0 GeV 



and jj duster i,2 < 1.6. Data and MC events are normalized in the region rj < 1.6 
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Figure 10: Energy (a) and p^ uster (b) distributions for third cluster: ZEUS data (full 
circles) and MC direct +resolved (full histogram). Data and MC are normalized to the 
number of events in data. 
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